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The Fermi liquid theory formulated by Landau is a basic paradigm of the behavior of an interacting 
many-body system. We present a new application of this theory to calculate "Landau force" on a 
macroscopic object. We show that immersing a pendulum in Fermi liquid can increase its oscillation 
frequency, and evidence of this has been observed in mixtures of 3 He and 4 He. 



PACS numbers: 67.10.Db, 62.10.+S, 67.60.G- 

Pendulum is one of the simplest but also most pre- 
cise physical instruments. For accurate measurements, 
one needs to understand the effect of the medium, often 
air, on the oscillation. The first effect perceived was to 
correct the gravitational restoring force by the buoyancy 
of the oscillating body. The second effect, calculated by 
Poisson in 1831, was that the flow of the medium around 
the body increases the inertia of the pendulum. The third 
correction, calculated by Stokes in 1850, was the effect of 
the viscosity of the medium pQ. With these corrections, 
the pendulum in fact turns out to be a device, a viscome- 
ter, to measure the properties of the medium. Further 
corrections are needed in a rarified gas, where the mean 
free path of the gas particles becomes comparable or ex- 
ceeds the size of the pendulum [2, 3 . In this article we 
continue this line of corrections by considering a pendu- 
lum in Fermi liquid. We show that the Fermi liquid acts 
like an elastic medium, which can increase the restoring 
force and thus lead to increased frequency of oscillation. 
We show that this effect is important in explaining os- 
cillator experiments made in liquid mixtures of 3 He and 
4 He [UH]. Moreover, this fermion-boson mixture allows 
a deep insight into the Fermi-liquid theory. 

The Fermi-liquid theory is a paradigm of what can 
happen in a strongly interacting many-particle system 
[BHS]. Instead of strongly interacting particles, the en- 
ergy spectrum has low-energy excitations called quasi- 
particles, which interact only weakly. Originally Landau 
formulated the Fermi-liquid theory for liquid 3 He [B], but 
the principle has much wider application, an important 
example being the conduction electrons in metals. The 
generalization of Fermi-liquid theory forms the basis to 
understand the superconductivity of many materials and 
the superfluidity in fermion systems. Even when it is not 
valid, it forms the standard against which to compare 
more sophisticated theories [9[ [10] . 

In the originai realm of 3 He, the Fermi-liquid theory 
was applied to explain many properties of the bulk, in 
particular thermodynamics, collective modes and trans- 
port properties [7[ [5]. Soon the theory was applied to 
a liquid limited by an oscillating planar wall [TT] IT^] . 
This problem stili had an analytic solution although it is 
very complicated, and therefore also approximate meth- 
ods are useful [12] ■ The approaches to non-planar geome- 



tries have been perturbative. As the Fermi-liquid theory 
reduces to the hydrodynamic theory in the limit of small 
mean free path, the leading corrections appear in the 
boundary conditions. This leads to a modified boundary 
condition where the fluid velocity is assumed to extrap- 
olate to the wall velocity at a "slip length" behind the 
wall [LllU8| . In this letter we calculate the response of 
the Fermi liquid to a vibrating body in the full range 
of mean free paths from the hydrodynamic to the ballis- 
tic limit. The calculation necessarily is numcric. With 
a view to apply the calculation to measurements using 
vibrating wires, we specifically consider a circular cylin- 
der oscillating transverse to its axis. We also use the 
Fermi-liquid theory adapted to the simultaneous pres- 
ence of bosonic superfluid, as formulated by Khalatnikov 
[15] . Besides being applicable to mixtures of 3 He and 
4 He, this has the advantage that some basic properties 
of the Fermi-liquid theory are more easily visible [20 . 

A noninteracting Fermi system has plane-wave states 
with momenta p and energy e p = p 2 /2m where m is the 
mass of a fermion. The ground state of the system con- 
sists of a Fermi sphere with ali states p < pp filled and 
others empty. The basic assumption of Landau was that 
the energy spectrum of the interacting Fermi system has 
one-to-one correspondence with the non-interacting sys- 
tem. In particular, the momenta p of the quasiparticles 
are the same as for noninteracting particles, but the en- 
ergies are shifted. Near the Fermi surface one writes the 
energy as linear in momentum, e p = vp(p — pp) where 
the parameter vp — pp j m* defines the effective mass 
m* . Using the group- velocity argument to the dispersion 
relation e p gives that the excitation propagates with ve- 
locity v whose magnitude equals the Fermi velocity vp. 
Thus the momentum of the principal fermion mv differs 
from the momentum p of the excitation. In Fermi-liquid 
theory this missing momentum is parameterized so that 
fraction D = 1 — (1 + |i*i)m/m* of the total momen- 
tum p is carried by the bosons and fraction mFx/3m* by 
other fermions. This can be interpreted that the princi- 
pal fermion pushes with it a cloud of both bosons and 
other fermions. Here the bosons are assumed to be fully 
condensed and are described by density pb, chemical po- 
tcntial fiB and superfluid velocity v s . The special case 
without bosons is obtained by setting p# = D = 0. 
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A cruciai observation of Landau was that the disper- 
sion relation assumed above leads to a consistent theory 
only if one allows an interaction between the quasipar- 
ticles. This modifies the quasiparticle energy to e p = 
vp(j> — Pf) + £>£p with the correction [TH] 

5ep(r, t) = (1 + a)5fi B (r, t) + Dp F p ■ v s (r, t) 

OC 

+ Y / Fi(Pi(P-P , )<f>p'(r,t))i } ,. (1) 

1=0 

Close to the Fermi surface, Sep depends on momentum 
only through its direction p = p/p. The first term in 
Eq. is the energy shift due to a non-equilibrium bo- 

son chemical potential <5/ìb = A'b — Mb^' parameterized 
by a. The second term is the energy shift due to mo- 
tion of the boson part with velocity v s . In the last term 
<f>p is the energy shift of the Fermi surface, Pq{x) = 1, 
P\{x) = x,. . . are the Legendre polynomials and (. . .)p is 
the average over the Fermi surface. The terms with coef- 
ficients Fq and Fi are analogous to the 1 + a and D terms 
but correspond to the fermion background. The higher 
order terms with l > 1 extend the leading two terms to 
general deformations of the Fermi surface. 

In order to determine the dynamics, one needs a kinetic 
equation. For small deviations from equilibrium it takes 
the form 



dt 



+ vfP ■ V( 



6e p ) 



I. 
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where / is the collision term. We see that this equation is 
a first order differential equation along classical particle 
trajectories, which are straight lines in the momentum di- 
rection p. The momentum conservation allows to deter- 
mine the stress tensor and thus the force on macroscopic 
objects. 

In order to see how the elasticity of the liquid arises, let 
us consider a beam of quasiparticles, as depicted in Fig. 
[T] On a trajectory crossing the beam, the beam causes a 
potential Sep, which shifts ali quasiparticle energies (Fig. 
[TJd). As in equilibrium the states are filled up to the Fermi 
level, the potential is compensated by a change of the 
particle density on the crossing trajectory. Therefore, the 
dynamics on the crossing trajectory, which is determined 
by quasiparticles at energies close to the Fermi level, is 
not essentially affected by a stationary beam. 

Let us now consider that the intensity of the beam is 
changing. The changing of the potential stores or re- 
leases quasiparticles on the crossing trajectory. Thus a 
varying quasiparticle beam radiates quasiparticles in ali 
directions even in the absence of any collisions (1 = 0). 
In the case of a pendulum, the body moving with velocity 
u generates in the radiai direction r a beam with ampli- 
tude proportional to u f. Typically the Fo term is the 
dominant interaction terni in (JlJ and thus the potential 
Sep oc F u ■ f. The quasiparticles radiated back on the 
body are proportional to Slp oc FqÙ ■ r, where the dot 
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FIG. 1: Illustration of quasiparticle dynamics on a trajectory 
crossing a beam of quasiparticles. (a) The quasiparticle tra- 
jectory (with parameter s) and the beam depicted in r space, 
(b) The quasiparticle dispersion relations e p at three loca- 
tions on the quasiparticle trajectory crossing the beam. For 
a stationary 5e p the quasiparticles travel at Constant energy 
and are not essentially affected by the beam. (e) A temporal 
variation of the potential (Sè p ^ 0) leads to energy shift of 
quasiparticles. 



denotes time derivative. This results in an extra "Lan- 
dau force" F ex FqÌl the body has to exert on the liquid. 
A force of the form F = Mù is well known in hydro- 
dynamic fìow, and M can be interpreted as the mass of 
the fluid that is dragged with the body. The diffcrcncc 
in the present case is that Fo can be negative. Rather 
than thinking of a negative mass, a simpler interpreta- 
tion is that the fluid has elasticity leading to an increased 
restoring force and increased oscillation frequency of the 
pendulum. 

In pure 3 He the parameter Fo is positive. The Landau 
force leads to increased effective mass of any objects, in- 
cluding ions aerogel [35] and deformations of the 

surface of the liquid [23], which effect has not been ex- 
plored yet. The case of negative Fo is realized in mixtures 
of 3 He and 4 He. Below we calculate in detail the case of 
a circular cylinder oscillating transverse to its axis. 

The system of equations (JTJ- <[2j) was solved numeri- 
cally for 4>p(r) on a grid surrounding an oscillating cylin- 
der using relaxation-time approximation for the collision 
term. The force F by which the cylinder drives the liq- 
uid is conveniently expressed by mechanical impedance 
Z = Z' + iZ" defined by F = Zu, where we assume 
time dependence exp(— iuit). The dissipative Z' and the 
reactive Z" are plotted in Fig. || as a function of the 
mean free path l. Two cases are shown by solid lines: a 
cylinder of radius a in a large chamber and in a slab of 
thickness 16a. We see that the confinement has strong 
effect cspccially on Z" . 

In the limit £ — > 0, the viscosity can be neglected and 
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FIG. 2: The force F = Zu exerted by an oscillating cylinder 
(radius a, velocity u) on the surrounding Fermi-Bose liquid 
as a function of the mean free path l, (a) The real part Z 1 
and (b) the imaginary part Z" of the mechanical impedance 
Z. The main results are shown by blue solid lines. They 
correspond to diffusive boundary conditions in a slab cham- 
ber, which has two piane walls at x — ±8a and the cylinder 
in the middle oscillating in the x direction. Other lines dif- 
fer from the blue ones in the following respects: (green solid 
lines) a large chamber, (dotted lines) hydrodynamic approxi- 
mation in a large chamber [T][15], (dashed lines) 50 % spec- 
ular scattering at the oscillating cylinder, and (dash-dotted 
lines) quasiparticle-absorbing chamber walls. Four different 
contributions to the ballistic limit of the blue curve are in- 
dicated (see text). (c) Plot of Z' vs. Z" showing the same 
results as above and experimental data (orange) from Ref. 
Uat 5.6% 3 He concentration. The parameters of the calcula- 
tion uja/vF = 0.017, F = -0.28 and Fi = 0.155 were fixed by 
independent measurements. The large chamber is a coaxial 
cylinder of radius b — 80a with absorbing walls. 



the liquid behaves like an ideal fluid. This corresponds to 
pure reactance, Zideai = —ina 2 pujG. Here p is the liquid 



density and G— 1 = 0(a/b) 2 is a small correction caused 
by a finite chamber dimension 6. In the experimental case 
4J Z[' Acal / an-pp-F — —0.50. With increasing l the viscos- 
ity of the liquid becomes important. This leads to in- 
creasing dissipation and the effective mass M = —Z"/ui 
grows because the increasing viscosity causes more liquid 
to be dragged by the oscillating cylinder. The hydrody- 
namic description becomes insufficient when l approachcs 
a. Compared to the hydrodynamic theory (dotted lines), 
the dissipation grows less rapidly and the fluid mass cou- 
pled to the oscillator starts to decrease. In the ballistic 
regime £ ^> a the impedance saturates to a value inde- 
pendent of l. 

The fact that Z" in the ballistic limit exceeds its ideal- 
fluid value in the limit l — > can be analyzed as follows. 
1) The fermion component of the fluid dccouplcs from the 
ideal fluid flow corresponding to contribution mnp to be 
subtractcd from p in ^idoai- (Here up the fermion number 
density and m the bare fermion mass.) 2) As discussed 
above, part of the bosons are bound to fermion quasipar- 
ticles, and therefore the density Dm*n F /{! + \Fi) has to 
be subtracted from p. These two contributions together 
constitute what is known as the "normal fluid density" 
p n = m*np/(l + §i*ì). 3) The Landau force, the main 
contribution coming from the negative Fq. 4) In a fi- 
nite geometry, there is an effect caused by quasiparticles 
reflected from the chamber wall back to the oscillating 
body. Ali these effects are shown separately in Fig. [2] 
We see that the four effects are ali in the same direction, 
and havc similar orders of magnitudo. 

Fig. [2]: gives a comparison of the theory to experimcnts 
[H [5] . Experimentally the mean free path was controlied 
by temperature, the lowest temperatures approaching the 
ballistic limit. We see that the confinement is cruciai 
in order to achieve agreement. The essential point in 
the data is that the frequency shift of the oscillator in 
the low temperature limit is positive and clearly larger 
than the normal fluid contribution (contributions 1 and 
2 together). This indicates the importance of the con- 
tributions of Fermi- liquid interactions and confinement, 
although their individuai contributions cannot be sepa- 
rated. 

In the calculation, the parameters of the slab chamber 
and the oscillating cylinder are fixed by the experiment. 
The Fermi-liquid parameters m* = 2.46, F — —0.28 
and Fi = 0.155 are based on measurements of specific 
heat and second sound velocity as analyzed by Corruc- 
cini [25]. The higher Landau parameters F[ with / > 1 
were assumed to vanish. The boundary condition corre- 
sponding to diffuse scattering of quasiparticles [TI] was 
used on ali surfaces for the blue solid lines in Fig. [2| Thus 
no fitting parameters was used. 

In more detailed comparison, one can allow for par- 
tial specularity of the quasiparticle scattering from the 
surface of the oscillating cylinder [TB]. Also, the walls of 
the experimental chamber are rough sinter, where part 



4 



of the quasiparticles likely is absorbed rather than scat- 
tered. The effect of these modifications are indicated by 
dashed and dash-dotted lines in Fig. [5] The effect of 
changing the other parameters and full comparison with 
the experiments [H [5] , which were made using two diffcr- 
ent wires and at different concentrations and pressures, 
will be presented elsewhere |26| . 

A Fermi liquid can support a collisionless propagat- 
ing mode, known as zero sound. The zero sound arises 
from the same physics as the Landau force discussed here 
[Eqs. Q-Q]. Using terms of antenna theory, the Lan- 
dau force could be interpreted as the near field effect and 
zero sound as the far field effect of the Fermi-liquid inter- 
actions. The zero sound occurs for repulsive interactions 
(Fq > in the case where F is the dominant interaction 
parameter). The elasticity effect (Landau force F and ti 
in opposite directions) occurs for attractive interactions 
(Fo < in the case of dominant Fq), where zero sound is 
damped. 

In the ballistic limit, the Landau force can be calcu- 
lated by low-frequency expansion of Eqs. ([l])-([2]). The 
result is 



7" 

Landau 



-a 2 m 



1=0 



QF, 



21 + 1 



Fi 



(3) 



The coefficients C'i ~ 1 are expressed in terms of compli- 
cated integrals. The dissipative part Z' in the ballistic 
limit has been studied in Refs. ITTI and IMI 

Before numerical calculation we made a transforma- 
tion that uncouples the Fermi and Bose parts in the low 
frequency regime |20| . which is a good approximation 
since ua/vp ~ 0.02. For this a new distibution function 

'tpp = 4>p + & e p + -Dto*^^-b/ to s(1 + \Fi) was defined. A 
two dimensionai grid was constructed around the cylin- 
der. On each grid point was stored the "old" values of 
the angular averages (ipp)p and (pippìp? which appear in 
and in / Then the kinetic equation was solved 
for ìpp(r) by integrating along the classical trajectories, 
taking into account boundary conditions [TTJ [5D]. For 
that the values of the angular averages had to be inter- 
polated between the grid points. By repeating this for a 
number of trajectories passing trough a given grid point, 
"new" values of the averages could be calculated by nu- 
merical angular integration. This was repeated for each 
grid point. Because of linearity of the equations, the 
new values of the angular averages are obtained as linear 
combination of the old ones. An inhomogeneous term ap- 
pears in this relation because of the boundary conditions 
at the oscillating cylinder. A self-consistent solution was 
obtained by requiring the new values to be the same as 
the old ones. This matrix equation was solved by numer- 
ical matrix inversion. 

In conclusion, we have shown that the force of a macro- 
scopic object on a Fermi liquid has a contribution from 
the interactions, caused in particular by F . With a nu- 



merical solution of the Fermi liquid equations in proper 
geometry we find good agreement with measurements in 
3 He- 4 He mixtures. 
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